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^— «J Using Molecular Dynamics simulations of a polymer liquid flowing past flat and patterned surfaces, we investigate the influence 
^| of corrugation, wettability and pressure on slippage and friction at the solid-liquid interface. For one-dimensional, rectangular 
! , grooves, we observe a gradual crossover between the Wenzel state, where the liquid fills the grooves, and the Cassie state, where 
O the corrugation supports the liquid and the grooves are filled with vapor. Using two independent flow set-ups, we characterize 
the near-surface flow by the slip length, 5, and the position, z^, at which viscous and frictional stresses are balanced according to 
^vq Navier's partial slip boundary condition. This hydrodynamic boundary position depends on the pressure inside the channel and 
£V| ma y be located above the corrugated surface. In the Cassie state, we observe that the edges of the corrugation contribute to the 
friction. 
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1 Introduction 

Understanding the flow of liquids past solid surfaces has at- 
tracted abiding attention in micro- and nanofluidics. While 
on macroscopic scales, the non-slip or stick boundary, which 
assumes the velocity of the fluid at the surface coincides with 
that of the surface, is popular, this approximation may become 
invalid on smaller length scales pertinent to the operation of 
micro- and nanofluidic devices^! Controlling the friction at 
the interface between the solid substrate and the liquid is cru- 
cial in these applications, and different strategies have been 
pursued: Like in the case of wetting, one can control friction 
by (i) the direct, microscopic interactions between the solid 
and the liquid or (ii) the surface topography of the solid 4 . 
While the former is largely dictated by the chemistry of the 
solid and liquid, the latter alternative is expected to be a uni- 
versal physical mechanism. 

We consider a regularly corrugated surface topography as 
sketched in Fig. [T] A liquid on this structured substrate may 
either fill the cavities between corrugations (Wenzel stated or 
build a straight liquid-vapor interface on the top of the corru- 
gations (Cassie state 6 or fakir state). The macroscopic contact 
angle that a liquid drop makes with the solid substrate is dic- 
tated by the balance of surface tensions according to Young's 
equation 7 . Surface roughness tends to amplify wettability, i.e., 
if the liquid did not wet the planar solid substrate, the contact 
angle would be even larger on a corrugated solid PEEI t 

The slippage of a liquid past a solid surface, in turn, is 
determined by a balance of viscous and frictional stresses at 
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Fig. 1 The geometries of the structured substrates. Patterning is 
present only in x-direction. 



solid-liquid boundary 10 . Modifying the topography of the sur- 
face, one simultaneously alters this balance and the surface 
tensions, therefore affecting both, slippage and wettability, re- 
spectively. A dramatical decrease of slippage takes place for 
a liquid in Wen zel state with respect to the corresponding flat 
substratd- 4 ! 11 ! 12 ! The effective position of equivalent no-slip 
wall was found between the top and bottom of the grooves^. 
In the case of the liquid in the Cassie (or fakir) state, in turn, 
the slip length increases and an enlargement by a factor of 2.5 
in comparison to flat substrates has been observed 4 . Semi- 
analytical models, assuming a no-slip condition at the solid- 
liquid interface and infinite slip at liquid-vapor interface, have 
been developed to calculate an effective slip! 12 ! 14 ! 

Both, in the description of wetting and slippage, the phe- 
nomenological approaches exploited the scale separation be- 
tween the geometry of the substrate (i.e., the shape of the cor- 
rugation) and the length scale that determines the surface ten- 
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sion or friction of the corresponding planar substrate. This 
assumption is justified for macroscopic corrugations but when 
the spatial scale of the corrugation decreases, geometry and 
interactions do no longer decouple; in this limit the surface 
roughness is an intrinsic property of the substrate. In this 
manuscript, we use Molecular Dynamics simulations of a 
particle-based model, which duly account for the discreteness 
of the fluid and its thermal fluctuations, to investigate the onset 
of deviations from the macroscopic phenomenological behav- 
ior. 

Keeping the fraction ip of the surface covered by grooves 
constant but varying their spatial dimensions, we investigate 
the influence of the roughness onto dynamical properties of 
the liquid in the Cassie and Wenzel states. Between these two 
states of the liquid, we only observe a gradual crossover but 
no phase transition 4 . We show that (i) the slippage at super- 
hydrophobic substrates is significantly influenced by the ad- 
ditional friction at the edges of corrugations and (ii) that the 
hydrodynamic boundary position does not necessary coincide 
with a localization of the solid-liquid interface between the top 
and the bottom of the grooves. 

Our paper is organized as follows: In section|2]we discuss 
the simulation model and techniques and the macroscopic phe- 
nomenological models. In this section we also provide infor- 
mation about the equilibrium wetting behavior. The compar- 
ison of channels with patterned and flat walls is reported in 
sections 3.1|3.3 The study of flows in channels with patterned 
walls at different pressures is presented in section 3.4 and we 
briefly summarize our results in section |4] 

2 Model and equilibrium wetting properties 

2.1 Coarse-grained polymer liquid and topographically 
structured substrate 

We represent the liquid by coarse-grained short chain 
molecules comprised of N ? = 10 monomersl 15 l 16 l The ex- 
cluded volume interaction and the cohesive attraction in the 
liquid are modeled by shifted Lennard- Jones (LJ) potential 
with energy scale e and bead diameter c. In the following, 
all energies are measured in units of e and all length scales 
are given in units of c. The LJ potential is truncated at r cut = 
2 x 2 1 / 6 er. Additionally, the neighboring monomers belonging 
to one polymer chain are bonded by FENE potentia l 17 ! 18 !. The 
virtue of using a polymer liquid is that the number density of 
the vapor phase is of the order of py ~ 10~ 5 c~ 3 , and one can 
neglect evaporation effects. Therefore the coexistence pres- 
sure is approximately zero. The density of the polymer liquid 
in coexistence with its vapor is p cocx = 0.786c~ 3 . 

The equations of motions are integrated by a velocity- 
Verlet algorithm using a time step of At = 0.005 in units 
of t — (jyjm/e, where m denotes the bead mass. A Dissi- 
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Fig. 2 (Color online) (a) Liquid in Cassie state does not fill the 
cavity, but creates straight liquid-vapor interface of area L y Wi. (b) 
Liquid in Wenzel state fills the cavity creating solid-liquid interfaces 
at the bottom and walls of the cavity. 



pative Particle Dynamics (DPD) thermostat) 19 l 2 ° l with friction 
coefficient 7dpd = 0.5 is used to control the temperature at 
k^T = 1.2e. The virtue of the DPD thermostat is the local 
conservation of momentum, which results in hydrodynamic 
behavior on large time and length scales. The self-diffusion 
coefficient D of the polymer liquid at coexistence pressure is 
D = 0.0157 ± 0.0030 cr 2 / 1 ^ Thus > the Rouse relaxation 
time of a polymer is tr = i? 2 e /(37r 2 £>) is tr = 25.3 ± 5 r. 

The flat solid substrate is represented by three layers of 
face-centered-orthorhombic unit cells of a LJ crystal with 
number density of p s = 2.67c -3 . The lengths of lattice vec- 
tors are a x = a y = \/3a and a z — a, where a = \/0.5a. The 
atoms of the substrate also interact with liquid particles by a 
truncated LJ potential II s (r) — Ul 3 (r) — U^(r c ) with 



(i) 



where the length scale is set by c s = 0.75er and the ratio e s /e 
is changed from 0.2 to 0.6 to control the wettability of the 
substrate. 

The topographically structured substrate has grooves as de- 
picted on a side view in Fig. [T] The ratio of the substrate cov- 
ered by grooves to the projected flat substrate is ip = 0.625. 
We study two structured surfaces with the same area fraction 
ip. To keep ip constant, we vary the distance between the 
grooves d and the width of the grooves w. These widths of 
the grooves are taken to be wi — 2.06<r and W2 — 6.18c. The 
height of the groove is identical for both substrates h = 2.38c. 

There are two microscopic length scales of the polymer liq- 
uid to which the dimensions of the substrate topography can 
be compared: (i) the end-to-end distance i?, ee = 3.43c of the 
spatially extended chain molecules and (ii) the effective bead 
size c. The former distance is related to the single-chain con- 
formations in a cavity (cf. Fig. [2|, while the latter sets the 
intrinsic scale of packing and layering effects in the fluid. For 
the short chain length considered in the present work, these 
two scales are not well separated. 
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The dimensions of the simulation box are L x — Ly — 33c, 
whereas L z is varied from 29.25er to 36.05cr depending of the 
topography of the substrates and strength of solid-liquid inter- 
action e s . If not mentioned otherwise, No = 1920 chains are 
confined in the simulation domain. The number of chains is 



varied in Sec 3.4 where we study flows at different pressures. 
The program code is written by authors. 

2.2 Wetting on the topographically structured substrate 

Structuring a hydrophobic substrate (contact angle of a drop 
greater than 90°) by grooves, one increases the contact angle 
and may render the substrate superhydrophobic. When the 
liquid does not fill the cavities (Cassie state), the contact angle 
changes according to Cassie's formula^ 



cos 8c — <p cos 6*e + f — 1 , 



(2) 



where 9c is the contact angle in the Cassie state, and #e is the 
contact angle on a flat substrate. 

In the Wenzel state, the liquid wets the substrate completely 
and the contact angle is described by Wenzel's formula^ 



cos 6v 



(3) 



where 6> w is the contact angle in the Wenzel state, and r stands 
for roughness of the substrate. The parameter r is the ratio of 
the actual substrate area wetted by the liquid to the projected 
area of the substrate, r > 1 and r = 1 for a flat unstruc- 
tured substrate. Our model of the substrate corresponds to the 
roughness r = d+ J^f h , where w is the width of the groove. 
For the finely and coarsely corrugated substrates the rough- 
ness is given by r\ = 1.865 and 7'2 = 1.289, respectively. The 
values of the contact angles for a flat surface (from anisotropy 
of the pressure tensor and Young's law^J) and the predictions 
of Cassie and Wenzel formulae (with roughnesses r% and r-i) 
are compiled in Table [T| 

We emphasize that Cassie and Wenzel predictions, given by 
Eqs. Q and Q, refer to macroscopic amounts of liquid and do 
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Table 1 Contact angles, 9e, 9c, and 6\y, on flat and structured 
substrates when the liquid is in the Cassie or the Wenzel states, 
respectively. Data are obtained from the measured values of the 
surface and interface tensions of the planar solid-liquid and 
liquid- vapor tension using Eqs. ^ and {3} for the two different 
corrugations r\ and r2 . 



not account for thermal fluctuations of the liquid-vapor inter- 
face spanning the grooves or the influence of the three-phase 
contact line at the edges of the grooves. Since the pinning of a 
liquid on a patterned substrate occurs in MD simulations, there 
is a hysteresis of the values of the contact angles. Thus, if one 
simulates a droplet, the estimates of the contact angle will not 
necessarily correspond to the Cassie or Wenzel values. 

To predict the stability region of the Wenzel state, we use a 
simple phenomenological continuum model^, accounting for 
free energies of the interfaces presented in Fig. [2] The liquid 
will fill the cavity if free energy of the liquid in the Wenzel 



state (Fig. 2b i is smaller than free energy of the liquid in the 



Cassie state (Fig. 2a I: 



F w - F c < 0. 



(4) 



From the geometries of interfaces the free energies can be 
written as 



F C = jLyWi 

F w = IShLyWi + j SL L y 2h = 



IShLyiwi + 2ft), 



(5) 



where Wi is the width of the cavity and 7sl is the solid-liquid 
surface tension that, in turn, is related to the contact angle 
of the liquid on a flat substrate 9% through Young's equation 
7sl = ~7 cos 6*e- Here we recall that the density of the vapor 
is negligible. Therefore we assume a substrate-vapor surface 
tension of 7sv ~ 0. Then, the condition of filling the cavity 
by the liquid can be formulated in terms of the contact angle 
6*e as 

9 E < arccos(— - — , — ). (6) 



1 + 2h/wi 

Substituting u>i with parameters of the substrates of our 
simulation model, we find that for a finely corrugated sub- 
strate the condition for filling the cavity corresponds to 6e 1 < 
108°, and for a roughly corrugated one it holds 6*e 2 < 124°. 
Comparing these phenomenological predictions with values 
of contact angles from Table. [T] we observe that the filling of 
cavities of finely and roughly corrugated substrates occurs at 
e s rj 0.6e and e s rj 0.5e, respectively. 

This phenomenological, macroscopic consideration pre- 
dicts a first-order surface phase transition between the Cassie 
and the Wenzel state. In Fig. [3] we present the number density 
of liquid in a single cavity of the substrate for the two topogra- 
phies. If there were a phase transition of first order between 
the Cassie and Wenzel states, one would expect a bimodal his- 
togram, where the two peaks corresponded to empty cavities 
without liquid (Cassie state) and cavities completely filled by 
the liquid (Wenzel state), respectively. In marked contrast, 
however, we observe only a single peak in the density distri- 
bution. As we vary the solid-liquid interactions, e s , the aver- 
age number density of liquid beads inside the cavity changes 
gradually. Around half-filling {p coex /2 ps 0.4ct~ 3 ) the distri- 
butions are very broad. The coexistence density of the liquid is 
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Fig. 3 (Color online) Histogram of a number density in a single 
cavity for finely (a) and roughly (b) corrugated substrates at varying 
strengths of solid-liquid interaction e s with body force 0.001a/r 2 . 
Additionally plotted a histogram in the crossover regime at 
e s = 0.4e without body force (dashed lines). The presence of a body 
force exerts only negligible influence on the density distribution. 
The coexistence density p COC x is plotted by dashed vertical line. 



indicated in the graph by the dashed vertical lines. The value 
inside the filled cavity is comparable but it remains smaller 
due to packing (layering) of the particle fluid inside the cav- 
ity. Snapshots reveal that the liquid-vapor interface strongly 
fluctuates in this crossover region and adopts states between 
straightly spanning the top of the grooves, as expected for the 
Cassie state, and touching the bottom of the groove, as being 
characteristic for the Wenzel state. In passing, we note that the 
application of the body force to the liquid does not change the 
behavior (cf. solid and dashed lines in Fig. [3] in presence and 
absence of a body force at e s = 0.4e, respectively). 

There are two effects that contribute to the absence of a true 
surface phase transition in our simulations: (i) If the coupling 
between neighboring cavities were negligible because the cav- 
ities were widely separated, the number of particles in each 
groove would be an independent quantity. Then, since each 
groove would be a quasi-one-dimensional system, there can- 
not be a true thermodynamic phase transition. Instead, we 



would expect to observe large but finite-sized domains of filled 
and unfilled portions along the groove. This rational is im- 
portant for widely spaced and extremely long grooves, and 
this mechanism will eliminate all thermodynamic singulari- 
ties at finite temperatures. While a systematic study of the 
filling/emptying behavior as a function of the spatial exten- 
sion along the groove is beyond the scope of the present study, 
it is unlikely that this effect causes the absence of a sharp ther- 
modynamic transition for the rather short grooves in our sim- 
ulations because we do not observe multiple domains along a 
groove. 

(ii) The absence of bimodal probability distributions in 
Fig.|3]rather indicates that even in the limit of infinitely long 
grooves, the transition between the Wenzel and the Cassie 
state is not of first-order but continuous. Thermal fluctuations 
give rise to fluctuations of the particle number in a cavity due 
to the finite compressibility of the liquid and, more impor- 
tantly, due to fluctuations of local position of the liquid-vapor 
interface inside the groove. The latter are the analog of capil- 
lary waves in this confined geometry. For the geometry used 
in our simulation study, the height of the grooves is on the 
same order of magnitude as the intrinsic width of the liquid- 
vapor interface, and already small excursions of the interface 
position result in a significant change of the filling fraction. 
Even if h/a is not very small but the grooves are widely sepa- 
rated to be considered independent, we anticipate large fluctu- 
ations because the macroscopic consideration, Eq. ([6|, asserts 
that the Wenzel and Cassie state have equal free energies for 
9e — > 90°. In this case, the liquid-vapor interface can move 
homogeneously like a rigid plane up and down in the groove. 
In fact, these simple macroscopic considerations are in agree- 
ment with a non-mean-field analysi s 23 " 25 !. 

If the liquid is confined into a channel, the pressure might 
deviate from the vapor pressure, at which the liquid and vapor 
of vanishingly small density coexist. The associated effects 
will be discussed in the following section. 

2.3 Normal pressure 

We study a polymer liquid that is confined between two ap- 
posing substrates. A flat unstructured substrate, a finely cor- 
rugated substrate with roughness ri and a substrate with a 
rough corrugation, r 2 , are considered. Additionally, we vary 
the strength of solid-liquid interaction, e s . To adjust the nor- 
mal pressure to its coexistence value, p cosx s» 0, we change 
the distance L z between the two apposing substrates and 
calculate the corresponding pressure inside the liquid com- 
posed of N — 1920 chains. The resulting width, L c ° ex , is 
therefore a function of solid-liquid interaction strength e s , the 
corrugation r,, and the number of chains No, i.e., L™ ex = 
Lf ex (<*, n, N ). 

We calculate the normal pressure tensor component p n l 26 l 27 l 
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dividing the systems into slabs perpendicular to the direction 
of the substrates and using the approach by Irving and Kirk- 



wood 1 ^ 1 . The pressure is divided into a configurational term, 
which is proportional to the number density p(k) in a slab k, 
and a virial term: 



Pn(k) = k B T(p{k)} 



1 



E 



(*) z\ 



■ cos (Fj 



(7) 

where V s \ denotes the volume of the slab, Xij , yij , Zij are the 
distances between interacting particles i and j in x,y,z di- 
rections, respectively, and rf 7 - = xf 7 + yf 7 - + zf 7 . The angle 
between the force F, 



r ij ■■ ,j i ',y 

and the distance vector r. 



. „ 11 iv> v^^x MJ is given by 

expression Fy, r^, and its cosine can take the values ±1. An- 
gular brackets (• • • ) denote averages in the canonical ensem- 
ble. The sum z2i<j runs over particles i and j if a portion of 
line connecting them is located inside the slab k. The interac- 
tions between the solid substrate and the liquid particles also 
contribute to the virial^. Since the solid-liquid interaction is 
of finite range, however, they do not contribute to the pressure 
in the center of the confined film. Mechanical stability asserts 
that the pressure normal to the surfaces does not depend on 
position and thus equals the pressure at the center of the film. 

The dependence of the normal pressure on the distance, L z , 
between the walls is presented in Fig. |4] At small L z , the den- 
sity and pressure at the center of the film monotonously de- 
crease with increasing L z . At constant number of particles, 
we observe that the higher the attractive strength of solid- 
liquid interaction, e s , is the smaller is L c ° ex , because of the 
excess of the compressible liquid at the attractive substrate. 
L z is defined as the maximal distance between the corrugated 
substrates, i.e., measured from the bottoms of the corrugation. 
L c ° ex is larger for a finely corrugated substrate with rough- 
ness r\ than for the coarse grooves (filled vs. open symbols in 
Fig. [4}, because the polymers explore the wider grooves bet- 
ter, and thereby they increase the liquid excess at the substrate 
and reduce the pressure at the center. 

At large wall separation, L z > 32.5a, we find deviations 
from the expected monotonous dependence on L z due to the 
formation of vapor bubbles under tensile stress, p < 0. These 
cavitation events nucleate more readily at strongly attractive 
surfaces e s ^ 0.5e and for corrugated substrates. In those 
cases, the simulations did not necessarily achieve equilibrium. 

3 Hydrodynamic boundary condition 

3.1 Navier's slip condition and computational techniques 

In order to study flow past these substrates, we move the two 
apposing substrates into opposite directions with a small ab- 
solute velocity i^ all = 0.075a jr (Couette flow). The average 
shear rate is defined as a ratio of the velocity of the wall w^ all 
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Fig. 4 (Color online) Pressure of the polymer liquid at the center of 
the film. Panel (a) shows the results for a flat surface, while panel 
(b) presents the results for the two corrugated substrates. Filled 
symbols represent the corrugation ri and empty symbols correspond 
to T2- Dashed horizontal lines indicate the coexistence density and 
pressure, respectively. 



to the width of the channel L z and, in our simulations, it is of 
the order of 7 ~ 5 • 10 -3 r _1 . The product of the shear rate 
and the Rouse relaxation time, t r , of the polymers defines the 
Weissenberg number, Wi = 7TR. Since Wi « 0.13 <C 1, we 
do not expect the weak flow to significantly perturb the molec- 
ular conformations and, indeed, density profiles obtained with 
and without flow are almost identical (see e.g. Fig. [3]). 

When the fluid flows past the substrate, the velocity profile 
deviates from the macroscopic behavior, and the velocity of 
the liquid at the surface may not equal that of the substrate, 
i.e., the liquid slips past the surface. These surface-induced 
deviations from the macroscopic behavior can be incorporated 
into a macroscopic continuum description via boundary con- 
ditions. A popular choice is Navier's hydrodynamic boundary 
conditionESl 



1-[T||5 



dv x 
az 



(8) 



where r\ denotes the viscosity of the liquid, v x is the x com- 
ponent of the velocity of the liquid parallel to the substrate. A 
is the friction coefficient, and the slip velocity, v s = v x (zh), 
represents the velocity of the liquid at the surface. If ap- 
plied to the detailed microscopic velocity profile, this equa- 
tion would quantify the balance between the viscous stress in 
the liquid and the friction stress of the liquid slipping past the 
surface. In order to serve as a boundary condition of a contin- 
uum description (i.e., the Navier-Stokes equation), however, 
changes of the liquid structure at the surface are ignored, and 
1) is interpreted as the shear viscosity of the liquid in the bulk, 



r? = 5.3 ±0.1o-7 



By the same token, v x {z) is inter- 



preted as the macroscopic velocity profile that obeys the con- 
tinuum description and that we can extract from simulations 
by extrapolating the linear profile of Couette flow towards the 
substrate, z — z^. 

Thus, Navier's hydrodynamic boundary condition, Eq. (jHJ, 
parameterizes the flow past the surface by two effective ma- 
terial constants of the surface^EUE! the hydrodynamic posi- 
tion, Zh, at which the boundary condition to the macroscopic 
continuum description is applied and the slip length, S = r]/X. 
The flow of simple liquids over flat substrates is often de- 
scribed by the no-slip boundary condition, S — 33 , and the 
hydrodynamic position, Zh, often coincides with the location 
of the sharp solid-liquid interface, as intuitivel y expec ted. Fi- 

~~ W at 
and chemi- 



nite slip has been observed in complex liquids " 



specific substrates: superhydrophobic 
cally patterned onesPSEl 

Here we study the flow of a polymer liquid past corrugated 
substrates. Therefore we expect a finite slip length, 5, and 
the effective location, Zh, at which the hydrodynamic bound- 
ary conditions is to be applied, is not obvious. Since Navier's 
boundary condition features two independent parameters, we 
use two flow profiles: Couette flow, generated by moving the 
surfaces, and Poiseuille-like flow, generated by applying an 
external body force on the liquid. The linear and parabolic 
velocity profiles expected from the macroscopic continuum 
model for Couette and Poiseuille flow, respectively, are illus- 
trated in Fig. [5] 

We chose the origin of the coordinate system at the cen- 
ter of the film. The macroscopic velocity profile is linear 
v x (z) = az + b for Couette flow, and the constants a and 
b are determined from the boundary conditions v x (0) = and 
v x (Zc) = — u^aii' wnere Zq < is the position where ex- 
trapolation of the linear, macroscopic velocity profile equals 
to the velocity of the bottom wall — v^ all as shown in Fig. 5a 





(a) 



(b) 



Fig. 5 (Color online) (a) Couette flow created in the system with 
roughly corrugated substrates. Two walls are moved into opposite 
directions with absolute velocities t^aii' The origin of the z-axis is 
taken in the middle of the channel (black). The linear velocity profile 
is sketched by the arrows of different size. The interpolated profile 
reaches the velocity of the wall at the coordinate Zc- (b) The flow 
created by application of a body force / (white) onto the particles of 
the liquid. The parabolic velocity profile (arrows) is interpolated to 
the velocity of the wall (v^aii = cr/r) and reaches it at Zp. 



The resulting velocity profile is therefore 



Zr 



(9) 



According to the macroscopic Navier-Stokes equation, the 
velocity profile of the flow generated by an external body force 
has a parabolic shape, v ? x (z) — — p ^ z 2 + Az + B, where 
Pcoex = 0.786cr~ 3 is the number density of the liquid and 7/ 
denotes its viscosity. A = by symmetry. When we express 
the constant B by the boundary condition v^(±Z P ) = 0, we 
obtain: 



2V 



(Z 2 - z 2 ). 



(10) 



Z P is the coordinate, where the extrapolated parabolic velocity 
profile of the liquid reaches the velocity of the wall as illus- 
trated in Fig. 5b Velocity profiles are fitted by Eqs. |9]and 10 
at the center of the film excluding the near-substrate region. 
Applying the hydrodynamic boundary condition, Eq. |[8), to 
the macroscopic profiles, Eqs. |9]l and ( fT0] > we relate the ex- 
trapolated positions Zq and Z P , which we extract from the 
simulations, to the slip length, 5 = \J Z^ — Zp, and boundary 
position, z h = Z c + 

To generate Poiseuille-like flow, we impose body forces 
/ of strength 0.001(t/t 2 , 0.003cr/r 2 , 0.005ct/t 2 and 
0.008ct/t 2 . An external force injects energy into the system, 
which is removed by the thermostat, keeping the temperature 
constant. In order to preserve a linear response of the liq- 
uid, only the smallest force was applied at small values of e s , 
whereas for e s > 0.4e all four forces were used. 
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Fig. 6 (Color online) Distance of the position, Zh, at which Navier's 
boundary condition[8]is to be applied, to the innermost substrate 
atom Zq P for channels with fine and rough corrugations (lines with 
filled circles and squares, respectively) and for channels with flat 
walls (line with open circles). Positive values, Zh — z'q P > 0, 
indicate that the effective boundary position is closer to the center of 
the channel than any particle of the wall. 



3.2 Position of the hydrodynamic boundary 

Using the simulation data obtained from Couette and 
Poiseuille-like flows at p coex rs 0, we can independently deter- 
mine the slip length, S, and the position of the hydrodynamic 
boundary, z^. In Fig. [6] we present the distance between the 
boundary position Zh and the position, Zq P , of the innermost 
substrate segment, which is located on the tops of ridges be- 
tween the grooves. In case of a flat substrate (black line with 
open circles), the position Zq V refers to the innermost atomic 
layer of the substrate. It is interesting to notice, that even in the 
case of flat walls and small e s , the hydrodynamic boundary is 
found somewhat inside the liquid, but not directly at the solid- 
liquid interface. For most parameters at structured substrates, 
the hydrodynamic boundary is located 2 — 4 segment diame- 
ters inside the liquid above the top of the ridges between the 
grooves. This observation for our molecular fluid even holds 
in the Wenzel state, where the liquid enters the grooves. 

The position of an effective no-slip plane, where the extrap- 
olation of the hydrodynamic velocity field, vanishes - Zc for 
Couette flow and Zp for Poiseuille flow - will be located at an 
intermediate coordinate between the peaks and valleys of the 
substrate topography in the Wenzel state if the slip length, 5, is 
larger than 2 — 4 segment diameters. This result is compatible 
with lattice Boltzmann simulations of a continuum fluid past 
statistically rough surfaces^]. We emphasize, however, that 
the no-slip position depends on the flow and, hence, cannot be 
used to parameterize a boundary condition, as shown in Fig. [7] 
Moreover, upon reaching the limit 9% 90° the slippage can 
be less than 2 segment diameters (cf. section 3.3 i, and the 



section 

corresponded no-slip position shifts into the liquid. 

In Fig. [8] we present the mass-weighted, two-dimensional 
velocity profiles for corrugated substrates at e s = 0.4e. Addi- 
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Fig. 7 (Color online) Velocity profiles of Couette flow (circles) and 
Poiseuille-like flow with body force / = 0.005<j/r 2 (squares) at 
e s = 0.4e and rough corrugation are shown. The case 
corresponds to the crossover regime, when liquid enters the cavity. 
Smooth thick solid lines are the fits to the velocity profiles, 
measured away from the substrate. The positions Zc and Zp, where 
extrapolated macroscopic profiles reach the velocity of the wall, 
v wM = -0.075cr/r and vl(-Z P ) = for Couette and 
Poiseuille-like flows, respectively, are indicated. These no-slip 
positions depend on the type of the flow and therefore cannot be 
used as a boundary condition. 



tionally, we demonstrate the velocity profiles v x (noisy solid 
lines) and their macroscopic fits (smooth solid lines) measured 
away from the substrate to neglect distortion by the friction 
at the boundary. The position of the hydrodynamic bound- 
ary condition is also indicated in the graph (horizontal dashed 
line). In both cases, it is above the corrugation. Qualita- 
tively, Zh describes the crossover between two behaviors: For 
larger distances from the substrate, the velocity is strictly par- 
allel to the x-axis as predicted by the Navier-Stokes equation 
for Poiseuille and Couette flow past a macroscopic substrate, 
while closer to the substrate, the velocity field acquires a spa- 
tially periodic, perpendicular z-component. In agreement with 
detailed hydrodynamic calculations of flow past a corrugated 
substrate, the amplitude of the perpendicular velocity compo- 
nent will decay faster away from the substrate if the lateral 
period is smaller 14 . Additionally, we note that the deviations 
of the velocity profiles from the macroscopic fits begin at the 
vicinity of hydrodynamic boundary position. 

Finally, we mention that, at e s = 0.4e (crossover regime), 
the average velocity of the liquid in a cavity of the finely cor- 
rugated substrate nearly vanishes, whereas in case of a rough 
corrugation the liquid actually flows in the cavity. This effect 
may influence the local shear stress and the effective boundary 
properties. 

3.3 Slip length and surface friction 

The second parameter of the hydrodynamic boundary condi- 
tion is the slip length, 5, or, equivalently, the friction coef- 
ficient A = rj/5. If the macroscopic continuum description 
of the velocity profile were accurate up to the position of the 
hydrodynamic boundary, Zh, and the liquid properties did not 
change in the vicinity of the substrate, then we could inter- 
pret the friction coefficient as the proportionality coefficient 
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Fig. 8 (Color online) Mass flux, pv (arrows), in the xz plane at e s = 0.4e and W\ = 2.06a in panel (a) and at e s = 0.4e and wi = 6.18<r in 
panel (b) under body force / = 0.005(t/t 2 (Poiseuille-like flows). Both cases, (a) and (b), correspond to the crossover regime, when liquid 
enters the cavity. In the accompanying plots velocity profiles (squares) and the extrapolated hydrodynamic profile (thick solid line) are 
shown.Smooth solid lines are the fits to the velocity profiles, measured away from the substrate. The position of the hydrodynamic boundary, 
Zh, is indicated by horizontal dashed lines. 



between the microscopic friction stress and the local velocity 
of the fluid flowing past the substrate. 

In this case, one can compute the friction coefficient via 
a Green-Kubo relation for the autocorrelation function of the 
tangential force, F x , between the liquid and the substrate^, 

R s =x = dt {Fx{t)Fx{0)) 

where A = L x L y is the area of the substrate. The analysis 
of Barrat and Bocquet 30 suggests that A be proportional to the 
square of the attraction between liquid and substrate, A ~ . 
Previous simulations for a Lennard-Jones polymer liquid on a 
Lennard-Jones solid corroborated this relation 36 , and we also 
observe this dependence on e s for our flat substrate as demon- 
strated in Fig. [9] Since the orientation of the unit cell of the 
crystalline substrate structure differs from our previous work, 
the friction coefficient, Afl at , that we observe is smaller than 
for the previous substrate. Such an orientation dependence of 
friction has also been observed in other studies'^. 

In this subsection we aim to adjust the pressure to its coexis- 
tence value, p C ocx ~ 0, such that the liquid- vapor interface in 
the Cassie state horizontally spans the grooves. Then, one ex- 
pects that the friction in the Cassie state is reduced, A = (fiX^t, 
compared to its value on a flat substrate, where tp is the frac- 
tion of the substrate area in contact with the liquid 4 . 

In the Wenzel state, the macroscopic expectation is that the 
friction, like the surface free energy, is proportional to the mi- 
croscopic contact area, r, between liquid and substrate. Thus 
a finely corrugated substrate, r*i = 1.865 is expected to give 
rise to more friction (with respect to the projected substrate 
area) than a coarsely corrugated one, r 2 = 1.289. 
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Fig. 9 (Color online) The slip length 8 for a liquid at flat (open 
circles) and patterned substrates of n and r-z corrugations (filled 
cirles and squares, correspondingly) at coexistence pressure. The 
inset shows the slip length 5 as a function of inverse squared 
interaction strength e s . 

Fig. [9] presents the slip length, S, as a function of the 
solid-liquid attraction. The qualitative behavior is similar for 
flat and corrugated substrates. For strong attraction, the slip 
lengths is microscopic but upon approaching e s — > (purely 
repulsive substrate), the contact angle approaches 180°(drying 
transition) and a lubricating thin vapor layer intervenes be- 
tween the substrate and liquid giving rise to large slippage. 

When plotting A ~ 1/6 vs. e^, one expects to find that 
the simulation data are compatible with two linear relations 
with different slopes. One linear relation corresponds to the 
Wenzel state for large e s , and another linear relation with a 
larger slope characterizes the Cassie regime, where the sub- 
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Fig. 10 (Color online) Friction coefficient A for different 
topographies of the substrate. The line with diamonds in the inset 
corresponds to the friction at the flat substrate multiplied by the area 
fraction ip — 0.625 covered by grooves. It is smaller than measured 
friction coefficients for corrugated substrates due to the friction at 
the edges of the grooves. The crossover regions, separating Cassie 
and Wenzel states, are shown by shaded areas for the different 
corrugations. 



strate has less contact with the liquid. In our simulation, 
however, the crossover gradually occurs in the wide range of 
e s /e G (0.275; 0.6) or e s /e G (0.25; 0.5) for the finely and 
coarsely corrugated substrate. Only at e s = 0.6e or e s ^ 0.5e, 
respectively, we observe the system consistently in the Wen- 
zel state. Therefore, we cannot identify the linear behavior 
that marks the Wenzel state. 

In the crossover region, the friction coefficients of both sub- 
strates do not significantly differ as shown in Fig. 
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cording to macroscopic considerations, we expected the finely 
corrugated substrate to generate more friction because of the 
larger contact area with the fluid. In our simulations, how- 
ever, the liquid is rather trapped inside the finely corrugated 
grooves, and the friction is not predominately generated at the 
solid-liquid interface but by viscous dissipation due to the ve- 
locity gradient at the interface between the trapped liquid in- 
side of the groove and the flowing liquid at the center. 

In the Cassie state, the friction of the corrugated substrate is 
comparable to that of the flat substrate. The macroscopic pre- 
diction that the friction is reduced by a factor of tp — 0.625 is 
not observed in our model. Moreover, the friction on the finely 
corrugated substrate is slightly higher than on the coarsely cor- 
rugated one. This effect can be partially rationalized by the 
effect of edges. 



In Fig. 1 1 we present two-dimensional density plots in the 



vicinity of the corrugated substrate. We observe that density 
oscillations that indicate layering effects of the particle fluid 
emerge from the edges of the grooves and propagate a few 
particle diameters, a, into the channel. We also note that even 
for the rather small body force, / = 0.005ct/t 2 , that drives 




Fig. 11 (Color online) Two-dimensional map of the number density 
of the particle liquid in the vicinity of the finely (left column) and 
coarsely (right column) corrugated substrate for e s = 0.3e (top row) 
and e s = 0.4e (bottom row). A body force / = 0.005(j/r 2 
generates a Poiseuille-like flow whose direction is indicated by the 
arrow. 



the fluid through the channel, the flow breaks the left-right 
symmetry, and the packing effects at the left and right edges 
of a groove differ slightly. The edge facing the flow (on the 
right side of the cavity) exhibits stronger layering effects than 
the edge on the left side. 

To a first approximation, we try to quantify these deviations 
from the macroscopic prediction due to localized perturba- 
tions of the flow by augmenting the Cassie model of friction 
by a term that accounts for additional friction at the edges of 
the grooves. 



A = <ysA flat + K(e s , d, w)S e d ge , 



(12) 



where S e d ge = is me li ne density of edges and 

K(e s ,d,w) quantifies the additional friction coefficient per 
edge length. K slightly depends on the geometry but in the 
limit of isolated edges, d,w — > oo, we expect K to adopt 
a finite value. The line density of edges is three-fold higher 
for the finely corrugated substrate than for the coarsely corru- 



gated one. From the data of Fig. 10 we extract the effective 
edge friction K. The values are in the range of 0.03 ^ K ^ 
0.16 for the finely corrugated substrate and somewhat larger, 
0.09 ^ K ^ 0.22, for the wider grooves. The order of mag- 
nitude of the edge friction, K ~ 2.5Afl a tcr, is plausible, i.e., 
each edge generates an additional friction that corresponds to 
roughly two additional rows of the substrate. The increase of 
K with increasing the distance between the edges is in accord 



with Fig. 11 where we observe that liquid layering effects 
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are more pronounced for wider grooves than for more narrow 
ones. 

3.4 Pressure-driven flow in patterned channels 

For an incompressible liquid confined into a channel with flat 
boundaries, the flow that is generated by a body force and 
pressure-driven Poiseuille flow are equivalent. For corrugated 
substrates, however, both methods of generating flow are no 
longer equivalent: If we use a body force to set the fluid in 
motion, the system will remain translationally invariant along 
the flow direction, and the pressure at the center of the chan- 
nel will be independent from the position along the channel. 
If we generate the flow by a pressure gradient, instead, the 
pressure will decrease along the channel. Even if the liquid 
were incompressible and hence its shear viscosity would not 
depend on pressure, we expect that the parameters of the hy- 
drodynamic boundary condition of a corrugated substrate - 
slip length, 5, and position, - will depend on the pressure 
and therefore will vary along the channel. The ability of the 
fluid to explore the cavities of the substrate makes the total 
system effectively compressible. Upstream, where the pres- 
sure is high, the fluid is more likely to enter the cavities of the 
substrate (Wenzel state), and the friction will be high. Farther 
downstream, in turn, the pressure is low and the liquid is more 
likely to adopt the Cassie state, which results in a lower fric- 
tion. At negative pressures, p < p coex , one might expect either 
again an increase of friction because the liquid- vapor interface 
in the Cassie state is curved (bubble mattress) 44 or a decrease 
of friction because the vapor layer at the solid is formed and 
acts as a perfect lubricant? 4 ^. 

In order to study how the slip length depends on pressure, 
we generate Poiseuille-like flow via a body force and vary the 
number of polymers inside the slit-pore. We measure the den- 
sity at the center of the channel and use the bulk equation of 
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state (EOS) to determine the pressure in Fig. 12 This pressure 
corresponds to the normal pressure inside the channel. 



The simulation results in Fig. 12 demonstrate that the slip 
length sensitively depend on pressure. If the pressure is high, 
the slip length, 5, is microscopically small and the friction co- 
efficient A is very high. In this case we observe a parabolic 
profile with an effective no-slip boundary condition. Upon 
reducing the pressure, the slip length increases. When ap- 
proaching the coexistence pressure, the slip length increases 
for weakly attractive substrates because the grooves are grad- 
ually emptied. At p coc ^ the liquid adopts the Cassie state. For 
stronger attractions, e s the liquid remains in the Wenzel state 
even when we approach p COC x- The continuous transition be- 
tween Cassie and Wenzel states at p COC x as a function of e s 
has been studied in Sec. |2.2| For p < p CO ex very large slip 
lengths can be achieved because a thin vapor layer is formed 
at the substrate that acts as lubricant. In this case, we observe 
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Fig. 12 (Color online) Slip length, S, as a function of the pressure 
measured from EOS of the channel confined by finely (a) or 
coarsely (b) corrugated substrates, respectively. The insets present 
the friction coefficient, A, and the vertical line marks p cocx . The 
widths of the channels L z at fixed e s and corrugation correspond 
to the widths L™ ex (e», n, No), calculated in Sec. [23 



plug flow because the slip lengths exceeds our system size, L z 
by far, and the friction is very small. Such a lubrication layer 
forms more readily for small e s . 

The dependence of the hydrodynamic boundary position, 



Zh, on pressure is shown in Fig. 13 As in the case of flows at 
coexistence pressure, discussed in Sec. 3.2 for most param- 
eters the hydrodynamic boundary is found inside the liquid 
above the top of the ridges between the grooves. The signifi- 
cant change of hydrodynamic position, zj, at low pressures and 
at finely corrugated substrate, which goes along with a large 
statistical uncertainty, may be explained by the formation of 
a lubricant layer of vapor at the solid boundary and the con- 
comitant, large increase of slippage (plug flow). 

Thus slippage and friction greatly depend on the pressure 
inside the channel, and the pressure dependence may signifi- 
cantly affect the flow profile across the channel and the pres- 
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Fig. 13 (Color online) Distance of the position, z n , at which 
Navier's boundary condition[8]is to be applied, to the innermost 
substrate atom z'^ p for channels with fine (a) and rough (b) 
corrugations, respectively. Positive values, z n — Zq P > 0, indicate 
that the effective boundary position is closer to the center of the 
channel than any particle of the wall. Different line types 
correspond to different number of polymers in a channel at fixed 
L z — L™ cx (e s , Ti, No), giving rise to different pressures. 



sure drop along the channel. Using Navier's hydrodynamic 
boundary condition, we compute how the flow rate, Q, de- 
pends on the slip length, 8, and half of the effective channel 
width, Zh 



= 2 



dy 



\{z)dz = - 2L ^ dp ^ 5{p) 



3rj dx 



1 



(13) 



where I 2 is the local drop of pressure along the channel. Mass 
conservation requires that Q be constant along the channel and 
the pressure profile along the channel is dictated by 



dp 
dx 



1 



1 



(14) 



If the slip length is constant along the channel, the pressure 
decreases linearly with the coordinate, x, along the channel. 
For the corrugated substrate, however, 5 increases with p and, 
consequentially, the pressure decreases rapidly upstream and 
more slowly downstream. In the case of negative pressures at 
the end of the channel (suction) and small e s , we expect plug 
flow with vanishingly small pressure drop. 

4 Discussion and Outlook 

In this paper we studied the flows of a polymer liquid in chan- 
nels with patterned walls by non-equilibrium Molecular Dy- 
namics simulations and compared the results to the channels 
with flat walls. First, we measured the normal component of 
the pressure in the bulk at different widths of the channels of 
different topography and compared the slippage and the fric- 
tion at the coexistence pressure. 



Significant differences between flat and corrugated sub- 
strates were found for liquids in Cassie state. Owing to the 
friction at the edges of the grooves, however, the friction co- 
efficient A does not scale like the ratio of area covered by the 
grooves, ip, to the projected area. For microscopically corru- 
gated substrates, a correction associated with the fiction at the 
edges of the corrugations is suggested, c.f. Eq. 12 Contrary 



to the expectation, the hydrodynamic properties of the liquid 
in Wenzel state are not significantly influenced by the sub- 
strate corrugation, r. For the parameters of our simulation, the 
difference with respect to a flat substrate is rather controlled 
by the area fraction ip, but not the size of the grooves. We 
show that there is no sharp phase transition between Cassie 
and Wenzel states but only a rather gradual crossover. These 
observations indicate that the macroscopic concepts cannot be 
straightforwardly extrapolated down to substrate topographies 
with dimensions that correspond to tens of fluid particles. 

Moreover, we illustrated that one configuration of the flow 
is insufficient for a measurement of the hydrodynamic bound- 
ary condition, i.e., a boundary condition to the macroscopic 
hydrodynamic continuum equations that describes the sub- 
strate properties independently from the type and strengths 
of the flow. Using Couette and Poiseuille flow, we extracted 
the hydrodynamic boundary position, Zb, and slip length, 8, 
that characterize the Navier slip condition. In applications to 
superhydrophobic substrates of complex topography the posi- 
tion of the hydrodynamic boundary is not intuitive. While the 
equivalent no-slip plane, which depends on the type of flow, 
might be located between the top and bottom of the grooves^, 
we find that the position of the hydrodynamic boundary, Zh, 
was almost always located above the top of the roughness. 

Finally, we modeled the propagation of the liquid through 
a long microfluidic channel with patterned walls accounting 
for the pressure dependence of the liquid morphology and the 
concomitant friction. We observed that by using a non-wetting 
liquid (contact angle on a flat substrate 9e > 150°) the rate of 
the frictional dissipation is significantly reduced on the whole 
length of a long channel in comparison to more wettable sub- 
strates. Additionally, we found that only for a liquid at low 
pressure (p < p co&x ) the microscopic corrugation affects the 
hydrodynamic properties, whereas for higher pressures it does 
not exert a pronounced influence onto the slippage and fric- 
tion. 
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